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In the framework of the fc t -factorization approach, the production of T mesons 
at the Fermilab Tevatron is considered, and the comparisions of calculated pr- 
distributions and spin alignment parameter a with the DO experimental data are 
shown. We argue that measuring the double cross section and the polarization of 
upsilonium states can serve as a crucial test discriminating two competing theoret- 
ical approaches to parton dynamics in QCD. 



1 Introduction 

Nowadays, the production of heavy quarkonium states at high energies is under in- 
tense theoretical and experimental study [TJ [SJ. The production mechanism involves 
the physics of both short and long distances, and so, appeals to both perturbative and 
nonperturbative methods of QCD. The creation of a heavy quark pair QQ proceeds via 
the photon-gluon or gluon-gluon fusion (respectively, in ep and pp collisions) referring to 
small distances of the order of 1/(2toq), while the formation of the colorless final state 
refers to longer distances of the order oil /[tuq a s (rriQ)]. These distances are longer than 
the distances typical for hard interaction but are yet shorter than the ones responsible 
for hadronization (or confinement). Consequently, the production of heavy quarkonium 
states is under control of perturbative QCD but, on the other hand, is succeeded by non- 
perturbative emission of soft gluons. This feature gives rise to two competing theoretical 
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approaches known in the literature as the color-singlet [3] and color-octet [I] models. Ac- 
cording to the color-singlet approach, the formation of a colorless final state takes place 
already at the level of the hard partonic subprocess (which includes the emission of hard 
gluons when necessary). In the color-octet model, also known as nonrelativistic QCD 
(NRQCD), the formation of a meson starts from a color-octet QQ pair and proceeds via 
the emission of soft nonperturbative gluons. The former model has a well defined appli- 
cability range and has already demonstrated its predictive power in describing the J/ip 
production at HERA, both in the collinear [5] and the £; t -factorization [5] approaches. As 
it was shown in the analysis of recent ZEUS [7] data, there is no need in the color-octet 
contribution, neither in the collinear nor in the fc t -factorization approach. 

The numerical estimates of CO contributions extracted from the analysis of Tevatron 
data are at odds with the HERA data, especially as far as the inelasticity parameter 
z = E^f,/E 1 is concerned [S]. In the fc t -factorization approach, the values of the color-octet 
contributions obtained as fits of the Tevatron data appear to be substantially smaller than 
the ones in the collinear scheme, or even can be neglected at all QUI EJ [12] . 

In the present note we want to compare the predictions of k^~ factorization ap- 
proach |13j with the DO experimental data on the pt~ distributions (14] and the polar- 
izaton of T (IS) meson states produced at the Tevatron energies [T5] . 

2 Numerical results 

In the fc t -factorization approach, the cross section of a physical process is calculated as 
a convolution of the partonic cross section a and the unintegrated parton distribustion 
T g {x, fcj., ii 2 ), which depend on both the longitudinal momentum fraction x and trans- 
verse momentum kx- 

o PP = J T g {x ll k\ Tl p?)T g {x 2 ,k\ T ,[i 1 )agg{xx,X2,^ (1) 

In accord with the fct-factorization prescriptions [El [17l [TH] US] , the off-shell gluon spin 
density matrix is taken in the form 

^=pX^/IM 2 = *^/IM 2 - (2) 

In all other respects, our calculations follow the standard Feynman rules. 

In order to estimate the degree of theoretical uncertainty connected with the choice 
of unintegrated gluon density, we use two different parametrizations, which are known to 
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show the largest difference with each other, namely, the ones proposed in Refs. [THl HH] 
and [20] . In the first case [16] , the unintegrated gluon density is derived from the ordinary 
(collinear) density G(x,/j, 2 ) by differentiating it with respect to [j 2 and setting /i 2 = fc|>. 
Here we use the LO GRV set [21] as the input colinear density. In the following, this 
will be referred to as dGRV parametrisation. The other unintegrated gluon density [20j 
is obtained as a solution of leading order BFKL equation [19] in the double-logarithm 
approximation. In the following, this will be referred to as JB parametrisation. 



D0 , Run 2 Preliminary, 1 .3 fb 1 




Figure 1: Predictions on the production of T(1S) mesons at the Tevatron: the dou- 
ble differential cross section (left panel), the spin alignment parameter (right panel) as 
function of pr- 



The production of T mesons in pp collisions can proceed via either direct gluon- 
gluon fusion or the production of P-wave states Xb followed by their radiative decays 
The direct mechanism corresponds to the partonic subprocess g + g — * T + g 
which includes the emission of an additional hard gluon in the final state. The production 
of P-wave mesons is given by g + g —> Xb only. All essential parameters were taken from 
our previuos paper [13]. Fig. 1 shows the comparison of our results for the T (IS) meson 
production with the DO experimental data (TU [15] . 

The calculations presented here are also valid for the T(3S) state, except the lower 
total cross section (by an approximate factor of 1 /3) because of the correspondingly lower 
value of the wave function | v I / t(3S)(0)| 2 = 0.13 GeV 3 . The state T (3S) is producted 
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by the purely direct production mechanism and the prediction on the spin alignment 
parameter a becomes less uncertain. 
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Figure 2: Predictions on the production of T(1S) mesons at the Tevatron: the differential 
cross section and a parameter as function of y. 



The polarization state of a vector meson is characterized by the spin alignment pa- 
rameter a which is defined as a function of any kinematic variable as 

a{V) = (da/dV - 3da L /dV)/(da/dV + da L /dV), (3) 

where a is the reaction cross section and <tl is the part of cross section corresponding to 
mesons with longitudinal polarization (zero helicity state). The limiting values a = 1 and 
a = — 1 refer to the totally transverse and totally longitudinal polarizations. We will be 
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interested in the behavior of a as a function of the T transverse momentum: V = |pr|. 
The experimental definition of a is based on measuring the angular distributions of the 
decay leptons 

dr(T-^+/O/dcos0 ~ 1 + acos 2 (9, (4) 

where 9 is the polar angle of the final state muon measured in the decaying meson rest 
frame. Fig. 1 (right panel) shows the comparison of our results for the T (IS) meson 
polarization with the DO experimental data [15] . The integration limits over rapidity 
were adjusted to the experimental acceptances of DO (|yr| < 0-6) at the Tevatron. The 
yellow band in Fig. 1 corresponds the NRQCD predictions, where the strong transverse 
polarization is connected with the gluon fragmentation mechanism. 

In Fig. 2 we show the rapidity distribution for T (IS) production in more wide rigion 
than the DO experimental data, and behaviour of the alignment parameter a as fuction 
of rapidity. We see that a becomes positive at large y. Therefore we propose to measure 
the double differential cross sections of quarkonium productions. 

In summary we have considered the production of T mesons in high energy pp colli- 
sions in the fct-factorization approach and compred the predictions on the spin alignment 
parameter a{pr) with new the DO experimental data. We point out that the purest 
probe is provided by the polarization of T(35) mesons. We proposed to measure the 
double cross section of quarkonium production and the spin aligment parameter in more 
wide kinematical region on the rapidity also. 
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